Effect of Nonmagnetic Impurity in Nearly Antiferromagnetic Fermi Liquid: 
Magnetic Correlations and Transport Phenomena 
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In nearly antiferromagnetic (AF) metals such as high-T c superconductors (HTSC's), a single 
nonmagnetic impurity frequently causes nontrivial widespread change of the electronic states. To 
elucidate this long-standing issue, we study a Hubbard model with a strong onsite impurity potential 
based on an improved fluctuation-exchange (FLEX) approximation, which we call the GV^-FLEX 
method. This model corresponds to the HTSC with dilute nonmagnetic impurity concentration. We 
find that (i) both local and staggered susceptibilities are strongly enhanced around the impurity. 
By this reason, (ii) the quasiparticle lifetime as well as the local density of states (DOS) are strongly 
suppressed in a wide area around the impurity (like a Swiss cheese hole), which causes the "huge 
residual resistivity" beyond the s-wave unitary scattering limit. We stress that the excess quasi- 
particle damping rate caused by impurities has strong k-dependence due to non-s-wave scatterings 
induced by many-body effects, so the structure of the "hot spot/cold spot" in the host system 
persists against impurity doping. This result could be examined by the ARPES measurements. 
In addition, (iii) only a few percent of impurities can causes a "Kondo-like" upturn of resistivity 
(dp/dT < 0) at low T when the system is very close to the AF quantum critical point (QCP). The 
results (i)-(iii) obtained in the present study, which cannot be derived by the simple FLEX approxi- 
mation, naturally explains the main impurity effects in HTSC's. We also discuss the impurity effect 
in heavy fermion systems and organic superconductors. 

PACS numbers: 72.10.-d,74.81.-g,74.72.-h,71.27.+a 



I. INTRODUCTION 

In strongly correlated electron systems, the presence of 
nonmagnetic impurities with low concentration can cause 
drastic changes of electronic properties of the bulk sys- 
tem. Thus, the impurity effect is a useful probe to inves- 
tigate the electronic states of the host system. In under- 
doped high-r c superconductors (HTSC's), nonmagnetic 
impurities (such as Zn) causes a huge residual resistiv- 
ity beyond the s-wave unitary scattering limit. More- 
over, NMR measurements reveal that both the local and 
the staggered spin susceptibilities are strongly enhanced 
around the impurity. Until now, the whole understanding 
of these impurity effect in HTSC have not be achieved in 
terms of the Fermi liquid theory. By this reason, nontriv- 
ial impurity effects in under-doped HTSC's are frequently 
considered as the evidence of the breakdown of the Fermi 
liquid state. However, similar impurity effects are ob- 
served in other strongly correlated metals such as heavy 
fermion (HF) systems or organic superconductors, near 
the magnetic quantum-critical-points (QCP). Therefore, 
we have to develop previous theories of impurity effect to 
see whether these experimental results can be explained 
in terms of the Fermi liquid theory or not. 

In HTSC's without impurities, various physical quan- 
tities in the normal state deviate from the conventional 
Fermi liquid behaviors in usual metals, which are called 
the non- Fermi liquid (NFL) behaviors. Famous exam- 
ples of the NFL behaviors would be the Curie- Weiss like 
behavior of \/T\T and the T-linear resistivity above the 
pseudo-gap temperature, T* ~ 200K. One of the most 
predominant candidates is the Fermi liquid state with 



strong antiferromagnetic (AF) fluctuations QIIEIEELOl- 
In fact, spin fluctuation theories like the SCR theory 
|2| and the fluctuation-exchange (FLEX) approximation 
IE IE IE 111 E3 have succeeded in explaining various 
NFL behaviors in a unified way. Recently low-energy ex- 
citations in connection with HTSC's were studied in ref. 
[TT|. These approximations satisfy the Mermin- Wagner 
theorem (see Appendix A). Moreover, pseudo-gap phe- 
nomena under T* are well reproduced by the FLEX+T- 
matrix approximation, where self-energy correction due 
to strong superconducting (SC) fluctuations are taken 
into account |lll5lll2l|. 

One of the most remarkable NFL behaviors in HTSC's 
would be the anomalous transport phenomena. For ex- 
ample, the Hall coefficient Rn is proportional to T -1 
above T* , and \R u\ ^ > 1/ne (n being the electron filling 
number) at low T |13| . Moreover, the magnetoresistance 
Ap/po is proportional to RuJpq oc T~ 4 , which is called 
the modified Kohler's rule They were frequently 

cited as strong objections against a simple Fermi liquid 
picture, because analyses based on the relaxation time 
approximation (RTA) do not work. However, recent the- 
oretical works have revealed that the current vertex cor- 
rections (CVC's), which are dropped in the RTA, play im- 
portant roles in HTSC's. Due to the CVC's, anomalous 
behaviors of i?n, Ap/p, thermoelectric power and Nernst 
coefficient are naturally explained in a unified way, based 
on the Fermi liquid theory [E HH 11(1 ITtI lisj . 

In the present paper, we study the effect of a single 
nonmagnetic impurity on a Fermi liquid with strong AF 
fluctuations. For that purpose, we developed a useful 
method of calculating the real space structure of the self- 
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energy and the susceptibility around the strong nonmag- 
netic impurity, on the basis of an improved FLEX ap- 
proximation. When the AF fluctuations are strong, we 
find that both local and staggered spin susceptibilities 
are enhanced around the impurity. Moreover, the resid- 
ual resistivity Ap per impurity, which is determined by 
the nearly parallel shift of p(T), can take a huge value 
beyond the s-wave unitary scattering limit. In addition, 
a "Kondo-like" insulating behavior (dp/dT < 0) emerges 
in the close vicinity of the AF-QCP. These drastic impu- 
rity effects in nearly AF systems come from the fact that 
the electronic states are modified in a wide range around 
the impurity, whose radius is approximately AF correla- 
tion length, £af- The present study provides a unified 
understanding for various experimental impurity effect 
in HTSC's, as universal phenomena in nearly AF Fermi 
liquids. Any exotic mechanisms (breakdown of the Fermi 
liquid state) need not to be assumed to explain them. 

We find that the simple FLEX approximation does not 
reproduce reliable electronic states around the impurity. 
For example, the spin susceptibility x s is reduced around 
the impurity. This failure comes from the fact that the 
feedback effect on the vertex correction for X s is neglected 
in the FLEX approximation. To overcome this difficulty, 
we propose a modified version of the FLEX approxima- 
tion, which we call the GV^-FLEX method. 

In the present work, we want to calculate the impurity 
effect at sufficient low temperatures, say 50K. For this 
purpose, we have to work on a large real-space cluster 
with a impurity site, which should be at least 64x64 to 
obtain the correct bulk electric state at lower tempera- 
tures. The FLEX approximation for such a large cluster 
is, however, almost impossible to perform numerically. 
Here, we invented the method to calculate the single im- 
purity effect in a large cluster, using the fact that the the 
range of the modification of electronic states due to the 
impurity is at least 5^6 except for in heavily under- 
doped systems. 



A. Previous Theoretical Studies 

Effect of a nonmagnetic impurity embedded in the 
Hubbard model or in the t-J model had been studied 
theoretically by various methods. The t-J cluster model 
(~ 20 sites) with a single nonmagnetic impurity was stud- 
ied using the exact diagonalization method H3- They 
found that the AF correlation is enhanced around the 
impurity. The same result is realized in two-dimensional 
Heisenberg models with a vacant site because quantum 
fluctuations are reduced around the impurity pH |22| . 
The same mechanism will account for the enhancement 
of the AF correlation around an impurity in t-J model. 
Moreover, in the t-J model, an effective long-range im- 
purity potential is induced due to the many-body effect. 
The authors discussed that non s-wave scattering given 
by the effective long-range potential could give rise to 
a huge residual resistivity per impurity, beyond the s- 



wave unitary scattering limit. It is noteworthy that an 
extended Gutzwiller approximation was applied for the 
impurity problem in the t-J model [23| . 

Also, the Hubbard model with a single nonmagnetic 
impurity was studied usin g th e random phase approxi- 
mation (RPA) in refs. [24l l25l l26| . By assuming an "ex- 
tended impurity potential" , they explained that the local 
susceptibility is enhanced in proportion to xq [oc T _1 ] of 
the host, reflecting the lack of translational invariance. 
Similar analysis based on a phenomenological AF fluctu- 
ation model was done [27|. The staggered susceptibili- 
ties is also enhanced due to the change of the local DOS 
(Friedel oscillation) [28f, However, the RPA could not 
explain the enhancement of local and staggered suscepti- 
bilities when the (5-functional) onsite impurity potential, 
which corresponds to Zn or Li substitution in HTSC's, 
is assumed. Therefore, results given by the RPA are not 
universal in that they are very sensitive to the strength 
of extended impurity potential. 

In the present work based on the GV^-FLEX approx- 
imation, we show that 5-functional onsite nonmagnetic 
impurity causes the enhancement of x s universally. As a 
result, onsite impurities cause a huge residual resistivity 
due to the nonlocal widespread change on the self-energy. 
We will see that the Gl^-method gives a unified under- 
standing of the impurity problem in HTSC. 



B. Experimental Results 

Here, we introduce several experimental results in 
HTSC's and the related systems, which we focus on in 
the present work. In later sections, we will discuss the 
origin of these experimental facts, and show that they are 
qualitatively well explained in a unified way as the effect 
of nonmagnetic impurities or residual disorders, on the 
basis of the GF / -method. 

(a) Magnetic properties: In optimally or under- 
doped HTSC's, a nonmagnetic impurity replacing a Cu 
site causes a localized momentum. A Curie like uniform 
spin susceptibility is induced by dilute doping of Zn in 
YBa 2 Cu 3 7 _ a: (YBCO). The Curie constant C per Zn in 
under-doped compounds (x ps 0.34) is much larger than 
that in optimally doped ones (x ~ 0) 29] . They also re- 
ports an interesting relation C cx Ap. Curie-like suscep- 
tibility was observed in Al-doped La2-,sSr,5Cu04 (LSCO) 
[30| . In Zn-doped YBCO compounds, site-selective 89 Y 
NMR measurements revealed that both the local spin 
susceptibility |.3lLI.'S2| and the staggered susceptibility |33J 
are prominently enhanced around the Zn-site, within the 
radius of the AF correlation length £af- The same re- 
sult was obtained by the 7 Li Knight shift measurement 
in Li-doped YBCO compounds _.34j, and by the 63 Cu 
NMR measurement in Zn-doped YBCO compounds [35| . 
These NMR studies show that the impurities does not 
trap holes, contrary to the suggestion by refs. [ilj l20|. 
This fact means that the impurity-induced local moments 
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result from the change of the magnetic properties of itin- 
erant electrons around the impurity sites. 

(b) Resistivity: Fukuzumi et al. observed p(T) 
in Zn-dopcd YBCO and LSCO for Zn concentration 
"imp = 0.02 ~ 0.04 In over-doped systems, the 

residual resistivity Ap(T) per impurity, which is deter- 
mined from nearly parallel shift of p(T) by impurities, 
is consistent with the value for 2D electron gas; p? mp = 
(4h/e 2 )n imp /n. However, Ap ~ (4S,/e 2 )n imp /|l - n\ > 
pP mp in under-doped systems. This fact suggests that 
the scattering cross section of Zn anomalously increases 
in under-doped HTSC's, as if an effective radius of im- 
purity potential grows due to many-body effect. In addi- 
tion, upturn of p(T) (dp/dT < 0) is observed below 50K 
in under-doped compounds. 

Such a prominent enhancement of residual resistivity 
Ap is also observed in HF compounds near the AF-QCP, 
which is realized under a critical pressure P c [37l |38| . 
Famous examples are CeCusAu (P c w 3.4GPa) |39| . 
CeRhIn 5 (P c « 2GPa) and CeCoIn 5 (P c ~ OGPa) 
pfl Ell ] . Note that the enhancement of Ap has nothing to 
do with the increase of the renormalization factor z (<S 1) 
by applied pressure, because Ap is independent of z in the 
Fermi liquid theory. In addition, the residual resistivity 
of an organic superconductor K-(BEDT-TTF) 4 Hg2.89Br 8 
(T c w 4K), which is close to the AF-QCP at ambient 
pressure, decreases to be about 10% of the original value 
by applying the pressure |4^. Such a drastic reduction 
of Ap cannot be attributed to the change of the DOS by 
pressure. 

Ando et al. have measured the resistivity in HTSC's 
under high magnetic field (~60T), which totally sup- 
presses the superconductivity 0, 0|. In LSCO, they 
found that the insulating behavior emerges under the 
original T c at B = 0. This upturn of p occurs when 
kpl ~ 13 1 in the a6-plane (I being the mean free 
path), so it has nothing to do with conventional local- 
ization in bad metals. Also, It will be independent of 
the opening of the pseudo-gap because the pseudo-gap 
temperature T* is much higher. Moreover, neither the 
weak localization or the Kondo effect due to magnetic 
impurities cannot be the origin of the upturn, because 
the (negative) magnetoresistance in the insulating region 
is independent of the field direction, and the insulating 
behavior persists under very high magnetic field. 

Sckitani et al. also found similar insulating behav- 
ior (dp/dT < 0) in under-doped electron-doped systems, 
M2-<5Ce,5Cu04 (M=Nd, Pr, and La), in the normal state 
under high magnetic field |45|. They discussed that the 
residual apical oxygens (about 1%), which works as im- 
purity scattering potentials, give rise to the insulating 
behavior. They expected that the Kondo effect occurs. 
However, the B-dependence of p(B,T) does not seem 
to be consistent with the Kondo effect. The upturn of 
p c along the c-axis is also observed in optimally doped 
Sm 2 -,$Ce ( 5Cu04 (8 w 0.14) in ref. 46]; the upturn is 
robust against the strong magnetic field (~ 45T), both 



for B || a and B || c. Recently, upturn of p in optimally 
doped PCCO had been observed |47j . 

(c) Local density of states: 

STM measurement |48| revealed that a single nonmag- 
netic impurity in optimally doped Bi2Sr2CaCu20s+5 
causes strong suppression of the superconducting state 
with radius ~ 15A. This "Swiss cheese structure" below 
T c had been suggested by the p-SR measurement [4J| as 
well as the specific heat measurement in Zn doped LSCO 
|50|. Reference reports that the radius of Swiss 
cheese hole increases as the carrier doping decreases. The 
radius is approximately £af, rather than the coherence 
length. This result suggests that the electronic properties 
in the Swiss cheese hole is strongly modified even above 
T c , because £af is a characteristic length scale in the 
normal state. Moreover, resent STM/STS measurements 
have revealed that the local density of states (DOS) in 
BiaSrzCaC^Os+.T; is very inhomogeneous in the atomic 
scale |5ll 152^ . The observed ununiformity originates from 
the weak scattering potentials from out-of-plane dopant 
atoms 51]. These experimental observations suggest that 
the DOS and the electronic states in HTSC's are quite 
sensitive to the disorder potential. 

II. FORMALISM 

In the present paper, we study a(NxN) square lattice 
Hubbard model with an impurity site: 

H = H + Pi m p, (1) 

H = ^ek4 CT c k CT + U^n^mi, (2) 

kcr i 

Pimp = I(noi + n oi), (3) 

where Po is the Hubbard model for the host system. In 
Pimp, I is the onsite nonmagnetic impurity potential at 
the origin (r = 0). Because the translational invariance 
is violated in the case of I ^ 0, the self-energy E(r, r'; e„) 
and the Green function G(r,r';e„) cannot be functions 
of r — r'. In the present paper, we concentrate on the 
strong impurity potential case (unitary scattering case); 
I — oo. We will take the / = oo-limit in the course of 
the calculation. 

We develop the method to study the impurity effect 
with strong potential on the basis of the FLEX approx- 
imation. When / 7^ 0, Green functions which compose 
the self-energy for / = get insertions of Ps in all pos- 
sible manners. Then, the self-energy is divided into two 
terms: 

£(!•.,, rj-;e n ) = S°(r l - vf,e n ) + 5E(rj, r 3 -; e n ), (4) 

where E° is the self-energy for 1 = 0, that is, the self- 
energy for the host system without impurity. <5£ repre- 
sents the cross terms between / and U . Here, 8T, does not 
contain any terms composed of only Ps. The full Green 
function for Pq + Pimp is composed of S°, 8T, and I. In 



4 



G 1 



I 



G 1 



SE 



G 1 



8E = E - E 



FIG. 1: Dyson equations for G 1 and G in the presence of I 
and <5E. 



the case of I = oo, ££(r,0) = -£E(0,r) = -£°(r,0) 
because E(r, 0) = E(0, r) = 0. The impurity potential 
also causes the nonlocal change in the self-energy, that 
is, <5S(r i ,r :) ) is finite even for Y i: Yj ^ 0. <5E(rj,i\j) will 
quickly converge to zero as or Vj are away from the 
origin. 

The Dyson equation for the real-space Green function 
in the matrix representation is given by 

G(e„) = G 00 (e„) + G 00 (e„)(E(e„)+i)G(e„) 

= G (e n )+G°(e n )(S±(e n )+i)G(e n ), (5) 

where = I&ifi&jfi represents the impurity poten- 

tial at the origin. G 00 (ri - ry; e n ) = S k (* e « + M ~ 
e k ) _1 e lk '( r4_I " is the non- interacting Green function, 
and G°(e„) = ([G 00 (e„)]^ 1 - S°(e„))- 1 is the interacting 
Green function without impurity (1 = 0). 
The Dyson equation is also written as 

G(e„) = G 7 (e n ) + G / (e„)^(e„)G(e„), (6) 
G\e n ) = G°( en )+G°( en )7G 7 (e„), (7) 

where G 7 represents the Green function composed of I 
and G°, that is, G 1 = G\st,=o- Equations © and @ are 
expressed in Fig^ 

In the present work, we calculate the self-energy for 
I based on the improved FLEX approximation, with- 
out taking any averaging with respect to the position of 
impurity. Hereafter, we propose the three versions of 
frameworks for calculating the self-energy as follows. 

(I) GF°-method : First, we introduce the "GI/°-FLEX 
approximation" , where the full Green function G is ob- 
tained self-consistently whereas any impurity effects on 
the effective interaction V° are neglected. Here, the self- 
energy is given by 

£[GV0](r»,r., ;e„) = T G[qvo] ( r « > r j > + e n) 
i 

x^(r„r j;Wl ), (8) 

where e n = (2n + 1)ttT and lui = 21 ■ ttT, respectively. 
Gjgvo] an d <55][gvo] = ^[gvo] — ^° satisfy the Dyson 
equation, ©. V° and E° are given by the FLEX approx- 
imation for the host system. Hereafter, we call eq. © 



the Gy°-method for simplicity, because the self-energy 
is symbolically written as G o V°. Here, V° and £° are 
given by 



y (r,,r j;W; ) = ^]TF°(q,^K q - (r *- r ^ 



(9) 



v°(q,w,) = £/ 2 Qx q >0 + ^x q >;)-n q (^) 

(10) 

x° s(c) (^) = KM • I 1 - , (ii) 

= -T^G° +k (^ + e„)G^ e „), (12) 

k,n 

S k (e„) = T^G k+q (e„+^ (q,^), (13) 

where G k (e„) = (ie n + M - e fe - £ k (e n ))- x . Xq s and Xq c 
are spin and charge susceptibilities, respectively, given by 
the FLEX approximation for the host system. 

Using the GU°-method, We can calculate the nonlocal 
change in the self-energy induced around the impurity, 
<5E. However, the nonlocal change in the spin suscepti- 
bility is not taken into account in the GV -method. To 
calculate this effect, we introduce other two methods as 
follows. 

(II) GU-method : Next, we explain the "GU-method", 
which is equal to the FLEX approximation in real space. 
In this method, both the self-energy and the effective 
interaction are obtained fully self-consistently. The self- 
energy in the Gl^-method is given by 



s [GV](ri,rj;e n ) = G [GV] (r i; r,-; uji + e„), 

l 

xU(r l ,r J -;cj ; ) 



(14) 



r2 /3„„, . 1 



vfa) = u z l^-r^i) + - n(^) 

(15) 

where Gjgv] an d <^[GV] — E[gv] — E° satisfy Dyson equa- 
tion © . The spin and charge susceptibilities in the GV- 
method, x s and \ c , are given by 

x s(c) = n(i-(+)tm)~\ (16) 

xG [GV] {rj,ri;e n ). (17) 

In the GF-method, the impurity effect on V is fully taken 
into account in terms of the FLEX approximation. How- 
ever, we find that the numerical results given by the GV- 
method are totally inconsistent with experimental facts. 
This is because the vertex corrections (VC's) for the spin 
susceptibility, which are dropped in the GU-method, be- 
comes significant in strongly correlated systems. We will 
discuss this point in §V\ 
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(III) GV^-method: To overcome the difficulty inher- 
ent in the GF-method, we propose the "Gy 7 -method" , 
where the Green function is obtained self-consistently 
whereas the impurity effect on the effective interaction 
is calculated in a partially self-consistent way. We will 
show that the Gl^-method is the most superior among 
(I)-(III). Here, the self-energy is given by 



^[GVI]( r i! r j; e n) 



T^G [GV i](ri,r 

I 

x^(r„r j;Ul ), (18) 

(19) 



where G[ GV i] and <5E[ GV rj 



^[Gvi] — S° satisfy Dyson 
equation (j^J. The spin and charge susceptibilities in the 
GV 7 -method, \ Is and x Ic , are given by 



= ft (l-i+jun 1 



(20) 

n^r^rj-;^) = -T^G 1 (y t ,vf,e n + uji) 

e n 

xG J (r i)ri ;e n ). (21) 

In x Is and \ Ic in the GV^-method, the self-energy cor- 
rection given by the cross terms between / and U [ST.] is 
not taken into account, whereas it is taken in the GV- 
method. Nonetheless, results given by the Gy 7 -method 
are completely different from results by the Gl^-method, 
and the former results are well consistent with experi- 
ments. For example, x Is given in ea. l|20|l is strongly en- 
hanced around the impurity, whereas x s given in eq. i|16|) 
is suppressed by ST. In fJ3 we will show that the lat- 
ter result is an artifact of the GF-method because the 
reduction of \ s due to ST is overestimated. As a result, 
the Gy / -method is much superior to the GV-method. 

Figure |21 expresses the self-energies for GV, GV 1 and 
GI^°-methods, respectively. These methods are equiva- 
lent to the FLEX approximation when 1 = 0. In later 
sections, we solve the single impurity problem in the pres- 
ence of the Coulomb interaction on the basis of these 
three methods. 



III. METHOD OF NUMERICAL 
CALCULATION 

In this section, we study the two-dimensional Hubbard 
model with an impurity potential at the origin. We work 
on the (N x iV)-square lattice with periodic boundary 
condition. N should be large (at least 64) enough to 
achieve the thermodynamic limit at low temperatures. 
On the other hand, the range of nonlocal change in elec- 
tronic states due to the impurity is only a few lattice spac- 
ings. Taking this fact into account, we calculate ST a ^ 
(a and /3 being the lattice points in real space) by the 



FLEX [1=0]: E° = 



.-. .-V 1 



G' 



-...-.V 1 



• GVO: E = 



A 

G 



• GV: £ 



A 

G 



V°= WG ] 



A A A 

V = V EX [G] 



• GVI: I 



S***^L A I A At 

= C J • V l = V FLEX [G ! ] 



A 

G 



A 

= V 



81=0 



FIG. 

netic 



5£ = £ -£ : cross term of U and I 

A An An / A A A 

G = G°+G°(I + 5L)G 



2: Self-energy for the Hubbard model with a nonmag 
impurity, given by the GV, GV 1 and GV°-method. 



(-N/2.N/2) 



(N/2.N/2) 



B 

(-M.M) (M,M) 

A 

impurity 
(-M-M) "(M-M) 



(-N/2-N/2) 



(N/2-N/2) 



FIG. 3: Schematic expression for region A and B. 



(improved) FLEX approximation only for |a|, \ f3\ < M 
(M < N), and we put SE ai p = for \a\ > M or \/3\ > M. 
Here we put M = 6 ~ 8, which is enough to obtain a reli- 
able numerical result. Here, we explain the method how 
to reduce the working area to a ((2M + 1) x (2M + 1))- 
square lattice in solving the single impurity problem in 
the (N x iV)-square lattice. This technique helps us to 
avoid considerable numerical difficulty. Hereafter, we use 
i, j, ■ ■■ (a, /3, • • •) to represent the lattice points in the re- 
gion A+B (region A) in Fig|3] 

According to eq.Q, G| • is given by 



Gi /^ii 



^0,0 



(22) 
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G° G° 

G{i = G° itj + I for i,j ? 0, (23) 

1 - iLr o,o 

where G°j(e„) = G°(r l ,r j ;e n ) = G°(ri ~ r j; e «) = 
j^Ek G k( e ") exp(ik-(ri-rj); G£(e„) is the Green func- 
tion given by the FLEX approximation without impurity 
potential. Therefore, G\ ^ in the limit of / = oo is 



Gqj — Gj .o — 0' 



3,0 



°"0,0 



(24) 

fori.j^O. (25) 



According to eq.©, the Dyson equation for the full 
Green function G a ^ inside of region A is given by 



G a ,p — G^g + A a ^,G 



0' 



P' 



(26) 



(27) 



Note that <5E Qj( 3 is finite only when both a and (3 are in- 
side of region A in the present approximation. By solving 
eo.lffify. we obtain 

G(e„) = (i-i(e n )) _1 G / (e„), (28) 

which is a ((2M + l) 2 x (2M + l) 2 )-matrix equation. In 
the limit of I — oo, G ail a is given by 



\ 



/3 ^<0 



(29) 



Note that G Qj o = Gq,/3 = 0. A a ^ in eq.(J25J) is given by 



A a ,p = D° at . 



ro Q,0> 

^0,0 



D a,f3 - ^G° 0,(5^,3-^. 

/3' 



(30) 
(31) 



Next, we study the spin and charge susceptibilities in 
the presence of an impurity. In the FLEX approximation, 
the equation for the magnetic susceptibility in real space 
is 



A+B 



iXi, 



(32) 



where i, j and I represent the lattice points in region 
A+B. Ilij is the irreducible susceptibility defined in 
ea. (|T7Jl in the GV- method or in eq.J2Tl) in the GV 1 - 
method. This equation is not easy to solve because of 
its huge matrix size; (A 2 x A 2 ). To solve this difficulty, 
we reduce the above equation to ((2M + 1) 2 x (2M + 1) 2 )- 
matrix equation (inside of region A) , approximating that 



EL.J — when i and/or j are in region B. Here, 

n? • is the irreducible susceptibility for the host system: 
n?_ 3 -(w,) = -TE e „ G? 7j (e n + w l )q_ i (e n ). 
Here, Xa (3 ms ide region A is rewritten as 



ia,0 = K^ + UY, 1 *-* 



7#0 



B 

i 

B 



(33) 



(34) 



The infinite series of summation in ea . (|34|l can be taken 
by using the following fact: In the case of / — 0, ea. l|33|) 
gives X ° a s -p = ^ Ekn°(l - ITI")- 1 exp(*k • (r Q - r/J )). 
As a result, we obtain that 



= n° i/3 + c/]Tn« i7X « 



7,0 



7 



(35) 
(36) 



where the summation of 7 contains the origin. By solving 
this ((2M + l) 2 x (2M + l) 2 )-matrix equation, <5x s is 
obtained as 



Sx = (x 0s - n° - cm V s ) (1 + ufT 1 • (37) 

Using eq. 137|l , the solution of eq. l-i-il) is given by 



X s = [1-UIL-USx ) (Tl + Sx j, (38) 

$Xa,f3 = 5 XaA l ~ S <x0fi), (39) 

where the factor 1 — 5 a /3,o m eq. (|39|l represents the elim- 
ination of 7 = in the summation in eq. I|33|) . Numer- 
ical calculation of eq. l|38[> is easy because its matrix size 
((2M + l) 2 x (2M + l) 2 ) is not large for M = 6 ~ 8. In 
the numerical study, we have to check that all the eigen- 
values of 1 — UIl — U Sx in ca. (|38|l are positive, because 
a single impurity in a paramagnetic bulk system could 
not induce a static magnetic order. 

In the same way, we derive the charge susceptibility, 
which is given by the following (A 2 x A 2 )-matrix equa- 
tion: 



1 ':) 



A+B 

I 



(40) 
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Taking the same procedure as used in deriving ea. (l3*51) . 
Xa a in the region A is given by the following ((2M + 
l) 2 x (2M + l) 2 )-matrix: 

i + un + usx ) (n + 6 X ). (4i) 

Here, <5Xa*/3 is given by 

tXa,p = 5x c a ,fs{^- 5 af 3fi), (42) 

4° = (x 0c - n° + ufi°x 0c ) (i - ufT 1 > (43) 

where = ^ E k n k (l + C/n")- 1 exp(zk • (r, - r,)) 
is the charge susceptibility for the host system. 

Finally, the spin and charge susceptibilities in the 
GI^ 7 -method, x Is and X /c , are obtained as follows: 

x ls = (i - uti 1 - us x s *) I (n / + 5 X S *), (44) 

x Ic = (l + UIl 1 + U5 X C *y 1 (fl 1 + §X C *) , (45) 

where II 7 is given in ea. H21|l . Note that we have to check 

that all the eigenvalues of 1 — UTI — USx hi cq. 1441) are 
positive in the numerical study. 

IV. NUMERICAL RESULTS FOR LOCAL DOS 
AND SPIN SUSCEPTIBILITIES 

In this section, we show several numerical results given 
by GV°, GV 1 and GV methods. In each methods, the 
self-energy is obtained self-consistently. We find that 
Gy 7 -method gives the most reliable results, irrespective 
that the fully self-consistent condition for the quasiparti- 
cle interaction is not imposed. In the GV-method, on the 
other hand, the reduction of % s due to ST [i.e., the non- 
local change in the self-energy given by the cross terms 
between / and U] is overestimated. In SjV] we will show 
that the reduction of V due to ST is almost recovered 
if one takes account of the VC due to the excess spin 
fluctuations induced by the impurity. By this reason, 
the GV^-method is the most reliable formalism. In the 
present section, we mainly focus on the numerical results 
given by the GV^-method. 

In the present numerical study, the dispersion of the 
conduction electron is given by 

£k = 2t(cos(k x ) + cos(ky)) + At' cos(k x ) cos(k y ) 

+2i"(cos(2fc x ) + cos(2fc y )), (46) 

where t, t' ', and t" are the nearest, the next nearest, 
and the third nearest neighbor hopping integrals, re- 
spectively. In the present study, we use the following 
set of parameters: (I) YBCO (hole-doped): to = — 1, 
h = 1/6, h = -1/5, U = 6 ~ 8. (II) NCCO (electron- 
doped): t = -1, t x = 1/6, t 2 = -1/5, U = 5.5. (Ill) 
LSCO (hole-doped): t = -1, t x = 1/10, t 2 = -1/10, 



U = 4 ~ 5. These hopping parameters are equal to those 
used in ref.|l5|. They were determined qualitatively by 
fitting to the Fermi surface (FS) given by ARPES mea- 
surements or the LDA band calculations. The shape 
of the FS's for (I)-(III), all of which are hole-like, are 
shown in ref. [T^. Because to ~ 4000K in real systems, 
T = 0.01 corresponds to ~ 40K. In the present numerical 
study, 64 x 64 k-meshes and 1024 Matsubara frequencies 
are used. 

The value of U used in the present study is rather 
smaller than the bandwitdh (IVband), irrespective that 
the real coulomb interaction is larger than Wband- This 
will be justified by considering that U used here is the 
effective Coulomb interaction U c s between quasiparticles 
with low energies. In fact, U c s ~ Wband according to the 
Kanamori theory based on the two-particle approxima- 
tion [53| . 

In Ref. |TH . the shape of the Fermi surface, the tem- 
perature and the momentum dependences of the spin 
susceptibility [x q (0)] and the quasiparticle damping rate 
[Im£q(— iS)] given by the FLEX approximation are ex- 
plained in detail. The obtained results are well consis- 
tent with experiments. For example, the q-dependence 
of x q (0) shows that £af = 2 - 3a in YBCO (n=0.85) 
at T = 0.02, which is consistent with neutron measure- 
ments. 



A. Local Density of States 

Figure 0] represents the density of states (DOS) for 
a hole-doped system (LSCO, n = 0.9) at T = 0.02, 
at the nearest-neighbor site (r = (1,0)) and the next- 
nearest- neighbor site (r = (1, 1)) of the impurity, respec- 
tively. Here, "host" represents the DOS without the im- 
purity given by the FLEX approximation, iVhost(w) = 
^Ek ImG k( w ~ iS )- 0n thc othcr h and, Nf iu) = 
—ImGi^u) — iS) at site I, where & is given in eq. 
In Ni(u>), effect of ST induced around the impurity is 
dropped. At (1,0) [at (1,1)], N[(lj) is larger [smaller] 
than A r host(w), which is recognized as the Friedel oscil- 
lation |28|. This result is changed only slightly by GV° 
or GV methods, which are not shown in However, 
the DOS given by the Gl^-method is much smaller than 
Ni (u>) in under-doped systems, because Im<5£ takes a 
large value around the impurity. The Green function is 
given by 

G [G vi] = ([G^-fetcvi])" 1 , (47) 

which are easily obtained by ea. ()29jl in the present nu- 
merical study. As shown in Fig. the local DOS given 
by the GT^-method is strongly suppressed, especially 
at (1,0). This suppression becomes more prominent for 
U = 5. The reason for this suppression is the extremely 
short quasiparticle lifetime, which is caused by the huge 
local spin susceptibility x Is around the impurity. We will 
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FIG. 4: (Color online) Obtained local DOS for LSCO (n = 
0.9, U = 5) at T = 0.02 (~ 80K). (i)host: DOS given by 
the FLEX approximation without impurity, (ii) N' (uj) — 

-^Ek ImG7 H where G ' 

is given in eg. 1251 . (hi) local 
DOS given by the GV 7 -method; —777 ImG[Gvi] (w) 
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NCCO (U=5.5. n=l'.15) 




I 2 3 4 5 6 



FIG. 5: The reduction of the local DOS at the Fermi level 
around the impurity given by the GV 7 -method, both for 
LSCO (n = 0.9, (7 = 4) and for NCCO (n = 1.15, U = 5.5) 
at T = 0.02, respectively. In the host, AW(0) = 0.322 for 
LSCO and AW(0) = 0.217 for NCCO, respectively. 



discuss the impurity effect on the spin susceptibility in 
the next subsection. 

Figure [3] shows the reduction of the local DOS at the 
Fermi level obtained by the GT^-method, both for hole- 
doped and electron-doped systems at T = 0.02. We see 
that the DOS is prominently suppressed in a wide re- 
gion around the impurity, especially along the diagonal 
axis. The suppression of the DOS is also recognized along 
x, y-axis, which is caused by the enhanced quasiparticle 
damping, Im<5£(— z<5), given by the GV^-method. The 
strong suppression of the DOS around the impurity site 
is consistent with the "Swiss cheese structure" observed 
in the STM measurement |48|. For LSCO, the radius 
of the Swiss cheese hole in Fig. is about 3a (a being 
the lattice spacing). It increases further for U = 5, ap- 
proximately in proportion to the AF correlation length 
£af(oc VT). Because quasiparticle lifetime is extremely 
short in the Swiss cheese hole, small number of impurities 
will induce the huge residual resistivity and the promi- 
nent reduction of T c in under-doped systems. In WIBI 



we will calculate the residual resistivity and confirm this 
expectation. 

B. Static Spin Susceptibilities 

Figure HJ1 shows the obtained local spin susceptibility 
X s (r,r) along (1,0) and (1,1) directions at T = 0.02. 
We see that x s ( r J r ) is significantly enhanced around the 
impurity as U increases, or as n approaches unity. The 
radius of area where x s is enhanced is about 3 ~ 4a, 
which would corresponds to the Swiss cheese hole in the 
DOS. This result is consistent with the impurity effect in 
under-doped HTSC's observed by NMR measurements. 
On the other hand, an opposite result is given by the GT^- 
method, which is not reliable as discussed above. Note 
that we have checked that all the eigenvalues of 1 — Uli — 
Uox in ea. l|38|l are positive in the present numerical 
study. 

Figure \7\ represents the uniform susceptibility 
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FIG. 6: (Color online) Obtained local spin susceptibility 
X s (r,r) for LSCO, YBCO and NCCO at T — 0.02. Around 
the impurity site, x s ( r i r ) i s strongly enhanced in the GV 1 - 
method, whereas it slightly decreases in the GV^-method. The 
former result is consistent with experiments. 



Uniform = N~ 2 EiJ Xft (P) for (n = 0.9, U = 5) 

given by the G^-method. ni mp is the concentration 
of the nonmagnetic impurity. Here, each impurity is 
assumed to be independent. We see that the uniform 
susceptibility without impurity decreases slightly at 
lower temperatures, which corresponds to the "weak 
pseudo-gap behavior" in HTSC's above the strong 




0.15 



FIG. 7: (Color online) T dependence of the uniform spin sus- 
ceptibility given by the GV^-method in the presence or the 
absence of impurities. O represents the Weiss temperature. 



pseudo-gap temperatures, T* ~ 200K. Surprisingly, 
Fig. shows that a nonmagnetic impurity induces an 
approximate Curie-Weiss like uniform susceptibility 
A X » n imp • fJ? cS /3(T +9). For U = 5, /i cff = 0.74 MB 
and = 0, which means that 42% of a magnetic 
moment of spin-i (fi e & — 1.73/xb) is induced by a single 
nonmagnetic impurity. A similar behavior is obtained 
for YBCO in the present study. The obtained induced 
moment is slightly smaller than the experimental value 
Meg ~ lfi B in YBa 2 Cu 3 6 . 66 (T c « 60K) In 
the present calculation, a simple relation A\ oc Xo 
predicted by previous theoretical studies [2J, l25l 1261 |27| 
approximately holds. The obtained /i c ff is the present 
study is much larger and consistent with experiments. 

Figure [S] shows the nonlocal spin susceptibility around 
the impurity given by the GT^-method, x s (r,r'). This 
result means that staggered susceptibility given by the 
GV^-method is strongly enhanced around the impurity. 
On the other hand, it is slightly suppressed in the GV- 
method. We consider that the former result is correct 
whereas the latter result is an artifact of the GV^-method 
because of the lack of the vertex corrections, (see flVl l 

Here, we discuss the origin of the enhancement of 
X /s : Reflecting the large N[(uj) at I = (1,0), the ab- 
solute value of n/ .(0) in the GV^-method is strongly 
enhanced especially for i = (1,0) and (a) j = (1,0), 
(b) j = (1,1), and (c) j = (2,1). In the case of LSCO 
(U = 4, n = 0.9) at T = 0.02, 11^(0) [n° (0)] becomes 
0.172 [0.160] for (a), -0.160 [-O'lOO] for (b), and 0.104 
[0.0784] for (c). This fact will give the enhancement of 
\ Is = II / (1 — Uli 1 )^ 1 around the impurity site since 
eigenvalues of (1 — Uli 1 ) become smaller. In contrast, 
an "on-site" nonmagnetic impurity does not give the en- 
hancement of spin susceptibility in the RPA |24l I2M |2(| . 
The reason would be that the strong reduction of IT ) 
due to thermal and quantum fluctuations are well de- 
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FIG. 8: (Color online) Nonlocal spin susceptibilities X s ( r i r ') 
for LSCO at T = 0.02 given by the GV and GV^-methods. 
Around the impurity site, the staggered susceptibility is en- 
hanced in the GV^-method, whereas it decreases in the GV- 
method. The former result is consistent with experiments. 

scribed in the FLEX approximation. Then, the reduc- 
tion of fluctuations due to an impurity would give rise to 
the enhancement of susceptibilities. 

In summary, we find that both local and staggered spin 
susceptibilities are increased around the impurity site, 
within the radius of about 3a (~ £af) at T = 0.02. Sim- 
ilar results were obtained for both the hole-doped sys- 
tems (YBCO and LSCO) and the electron-doped ones 
(NCCO), rather insensitive to model parameters. More- 
over, similar impurity effects were obtained in the t-J 
model, by exact diagonalization studies 0, |2(j. As a 
result, the obtained impurity effects in this section will 
be universal in systems close to the AF-QCP. 

V. FEEDBACK EFFECT AND VERTEX 
CORRECTIONS 

In previous sections, we show that the local spin sus- 
ceptibility given by the GV^-method, % /s , is promi- 
nently enhanced around the impurity. However, \ Is 
could be modified if we go beyond the G^-method be- 



cause the induced susceptibility around the impurity, 
A% = \ Is ~ X° s i changes the susceptibility itself, in the 
form of the self-energy correction and the vertex 

correction (VC). We call this self-interaction effect "the 
feedback effect" . In the feedback effect, the susceptibility 
is enhanced by the VC whereas it is reduced by <5£. In 
the GF-method, where only the latter effect is taken into 
account, the local susceptibility becomes smaller than the 
host's value. This inconsistent result suggests the impor- 
tance of the VC's. In the present section, we study the 
VC's for the spin susceptibility, and show that the VC's 
almost cancel the self-energy correction. We find that the 
total feedback effect is very small in the GV^-method with 

gvl u ?°= °<^>° 




FIG. 9: Irreducible susceptibilities given by the GV 1 , GV and 
GV+VC methods, respectively. 

VC's (Gy+VC-method) : The obtained spin susceptibil- 
ity is similar to that by the Gl^-method. Therefore, we 
conclude that the GT^-method is superior to the GV- 
method. 

The irreducible susceptibility given by Gl^-method 
(n / - <T<T ) and GF-method (Xl aa ) are expressed in eas. l(T7j) 
and (|2T|) . respectively. Hereafter, we discuss the feed- 
back effect for the irreducible susceptibility perturba- 
tively with respect to AV = V 1 — V°. In this respect, 
W a can be expanded with respect to AV as 



II?J(0) « -2T Yl G^^Gl^GlAen) ■ AV^ r (e n - e m )Gi,^,(e m ), (48) 

i',j',e m ,e„ 

up to the lowest order, which is expressed in Fig. El We have checked numerically that eq. (|48|l is satisfied well. 

Next, we study the VC for irreducible susceptibility up to the second-order with respect to AV. The lowest order 
term, which we call the Maki-Thompson (MT) term customarily, is given by 

U 2 T 2 

AnjJ T (i,i) = — F(i ) j,i' ) j';e m) e n )(Axi l>r {e m -e n ) + Axl >r (em-e n )), (49) 

i'j';e m ,e„ 
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Anft T (i,i) = - 



U 2 T 2 



FihJii' >j'\ e m,en)AXi>,j'{£m - e n ), 



i 1 ,j';e m ,e n 



Ax 



Is.c Os.c 



F(i, j, i', f; e m , £n) = Gj, !i (e m )Gf )i , (e m )Gj, j (e n )G?jj, (e„), 
which is shown in Fig. |5J We also discuss the second-order term given as 

An A £(u) =t £ F'C^H.ia^O^'^ji^^^O+^'aji,^;-^)) 

x (^w (« 1 1 *2 , ji , J2 ; ^ ) - X° a , (i i , i 2 , ji , ]2 ; uji) ) , 
F'(i,i 1 ,i 2 ;u l ) = T^Gf 2!l (e„)G^ (e„)Gf l!l2 (e„+^), 

X^(ii,i 2 ,ji ) j 2 ;^i) = —XiljMlh + ^Mlji( v X&to + —x\l n {Ux^ un +45 ilA ) 

7 ^iiij'i^»2iia ' 



(50) 

(51) 
(52) 



xl l (ii,i 2 ,ji,j 2 ;ui) = —(xf", n - xl^Hxt^ " *tjJ + ^Y s h,h(x! s 2tj2 - xll 32 ) 



u 3 



IP 



where £ = or /. We call ea. (|53|l the Aslamazov-Larkin term. 



(53) 
(54) 



(55) 



(56) 



In the FLEX approximation, the irreducible VC's given 
by the Ward identity, r lIT = 6T,/SG, are composed of the 
MT-term and the AL-term @|. 

As a result, the irreducible susceptibility given by the 
"GV- method with VC's up to the second-order with re- 
spect to Ax s ' c " (GV+VC method) is given by 

n™'(i,j) = n(i,j)^; 

+AIl™;(i,j) + AW A i(i,j). (57) 

The spin susceptibility in the GV+VC method is ob- 
tained as 



(3,0) 



Xgv+vc 

xtt 
xn 



xn+xu> (58) 
ti^ +tL^Uxn-n n Uxi h (59) 



In the present numerical study, we calculate 
AIImt,al(«, (3) only for region A in Fig. |3] Here we 
put M = 3 in the present section. Unfortunately, it is 
not easy to calculate all the elements of eqs. (|49f) . (|50f) 
and (|53|l in the region A because of the huge computa- 
tion time. Therefore, we calculate AIImt,al(Q!, P) only 
for | a — j3\ < 4 in region A, and derive the static spin sus- 
ceptibility Xgv+vc by solving ecis. l(5*8|) - i|6TJ|) for region 
A+B. 

Figure ^] show the local spin susceptibility given 
by GV 7 -method [x /s («,?)], GV-method [x s (i,z)], and 
GV+VC-method [xbv+vcl^ 0] for LSC0 (f7 = 4, n = 
0.9), respectively. At i = (1,0), we see that x /s (M) 
increases whereas X s (i,i) decreases, as we have shown 



(1,0)-(3,0) 
t— • GVI 
I BGV 

GV+VC (MT only) 
GV+VC (MT+AL) 




LSCO (n=0.9) 
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FIG. 10: (Color online) Local spin susceptibilities for U — 4 
at T = 0.02 given by the GV, GV 1 and GV+VC methods, 
respectively. V( r > r ) at r = (1,0) takes an enhanced value 
for both GV 1 and GV+VC methods, whereas it decreases in 
the GV-method. We conclude that GV / -method is as good 
as GV+VC-method. 



in Fig. El Note that we have put II(i,j) = Uj(i,j) = 
IIgv+vc^, j) = for \i — j\ > 4 in deriving Fig. HOI to 
make comparison between different methods. Because of 
this fact, both x /s (z, i) and x s (i, i) in Fig. Hj]]are smaller 
than those in Fig. [OJfor U = 4. 

As shown in Fig. Xgv+vc (M) at i = (1,0) is 
strongly enhanced due to the VC's, to be comparable 
with % /s (i,i). This enhancement is brought mainly by 
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the MT-term, whereas the AL-term slightly reduces the 
spin susceptibility. Therefore, the suppression of \ s in 
the GV-method, which is caused by the non-local self- 
energy correction <5£, is almost recovered by the VC's. 
Therefore, the GV'-method gives a reliable spin suscep- 
tibility around the impurity, because the total feedback 
effects almost cancel. 

We comment that the superiority of the GV-method 
over the GV-method and the importance of the VC's 
would remind us of the GW approximation [H3 . l55t |5(| . 
It is a first principle calculation for the self-energy, which 
is given by the convolution of the Green function G and 
the screened interaction W within the RPA. In the "fully 
self-consistent GW" , both G and W are obtained self- 
consistently. In the GWo scheme, on the other hand, 
the self-energy is given by G and Wq, where Wo is 
the screened interaction without self-energy correction. 
Reference 0, EE Ell shows that the descriptions of 
the bandwidth reduction and the satellite structure in 
the quasiparticle spectrum are satisfactory in the GWo, 
whereas results given by GW are much worse. This result 
clearly indicates the necessity of VC's in the GW. 



VI. TRANSPORT PHENOMENA 



I I G 1 

t = X + >v-^— X 



+ Xr 



m — x 

5E 



-d) 
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FIG. 11: Diagrammatic expression for the t-matrix given in 
eqs. (CD-GEl- 



V,fJ, 



(64) 

(65) 



where t^p ~ 1S schematically expressed in Fig. ^2 
a and (3 represent sites in region A. 

Here, we derive the expression for the t-matrix in the 
limit I — ► 00. First, the second term of eq. I|62[l is 
rewritten as 



In previous sections, we showed that the magnetic 
susceptibility is strongly enhanced around the impurity. 
This fact will cause the strong nonlocal change in the self- 
energy around the impurity, <5£. In the present section, 
we will show that i5S gives rise to a huge residual resis- 
tivity at finite temperatures in nearly AF Fermi liquids, 
which is the most important finding in the present paper. 
Moreover, a small number of nonmagnetic impurities can 
cause a "Kondo-like" insulating behavior (dp/dT < 0) at 
low temperatures, when the system is very close to AF- 
QCP. Different from a conventional single-channel Kondo 
effect, the residual resistivity can be much larger than 
the value for s-wave unitary scattering. These findings 
naturally explain various long-standing problems on the 
transport phenomena in HTSC's, heavy fermion systems 
and organic superconductors. 



A. T-matrix 

First, we derive the expression for the t-matrix, t(e), 
which is defined as G = G° + G°tG . Therefore, the t- 
matrix due to the impurity potential / at center r = (0, 0) 
is given by 



I 2 Gq.q — I 2 Gq q + I 2 Gp aSY^a^Gg p 



+I 2 E G iJ^ G ^^GU- (66) 



(67) 
(68) 



According to eq. 

i 2 Gl = -i--^ r + o(r 1 ), 
^0,0 

igL = - Jf^ + oii- 1 ). 

Thus, eq. (|62[l in the limit / —> 00 is given by 
1+1 G ,o - --^Q— + 2^ 77% \ 2 



x hE a , + J2^a,' r G^sS'Ss,0j (69) 
In the same way, eq. I|63l) in the limit I — > 00 is given 



by 



n=l 

*a fl = (I + I 2 Gofi)S a fiS[lfi, 



(61) 
(62) 



ta fi = 1^2, (Go,v8Eu,P ■ Su,0 + 8T, a<v G v fi ■ 6/3,0) (63) 



y~] n o ' S ( G° s .o + y~] Gs^SE^G^fi ] So,, 



in 



+ (a^ (3). 



(70) 



Note that both eqs. I|64|l and l|65|l contain / only through 
<5£, which does not diverge even when / = oo. 
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To study the effect of the impurity on the transport 
phenomena, we take the average of the t-matrix with 
respect to the position of the impurity. The obtained 
result is 

A 

Tl = ^""^ ti+u,a 
ot 

^T- ( 1 - 7W~ E D la (<W + D a ^)G% ) 5 l>0 
^0,0 ^ ^0,0 afj J 

(71) 
(72) 



-grEl (5 a , + D a ,) 



0,0 
A 



A 

' E/ ^l+a,a, 



7, a 



(73) 



(74) 



where eqs. 1(72 ^1 . l(7l )l and J73J come from eqs. (E2J), 

j53t ) lEH and jnSl), respectively. D° ^ is given in eq. 
(|3*T)) . Similarly, -D Q „3 is defined as 

A A 
Da, (3 = E G' cti7 (5I] 7i /3 = E <^/3,7^7,a- (75) 

7 7 

Note that Z? Q:0 = Z? ,a = and Q — 0, whereas 
£>0a 7^ 0. After the analytic continuation of T}(e„), the 
quasiparticle damping rate (without the renormalization 
factor) due to the impurity is given by [EH l58| 



imp/ \ 

7 k ( £ ) 



iV 



^ImT,(e-i,5)e' k ' r ', (76) 



where ni mp is the density of impurities. 

Figure ^| shows 7™ p (0) along the Fermi surface for 
LSCO at T = 0.02, obtained by GV°, GV and GV 1 - 
methods. Here, we put ni mp = 1- 7 £ ost = ImE^— iS) 
in the host system given by the FLEX approximation. 
Filled square and triangular represents 



7i° mp ( e ) 



Im 



''imp 



G° (e-i5Y 



for n imp = 1, where Gg(e) = J2k ( e ) * s tne l° ca l 
Green function of the host. Equation l|77|l is a well-known 
expression for the quasiparticle damping rate due to s- 
wave unitary impurities. In fact, 7;° mp is derived from eq. 
GHl by putting <5£ Q:(3 = • 5 af3fi in eas. l(7T|) - (f7ljl . 

When the particle-hole symmetry is approximately sat- 
isfied around the Fermi level, ea. ljTTjl becomes 7i° mp (e) = 
n imp /TrN host (e), where N host (e) = ImGg(e -iS)/ir is the 
DOS of the host system. 

In each method (GV, GV 1 and GV°), 7 ™ p has a 
strong k-dependence similar to 7 £ ost , as shown in Fig. IT21 




(71,0) 



(JI/2.JI/2) 



(0,71) 



FIG. 12: (Color online) k-dependence of 7j c mp (0) per ni mp at 
T = 0.02, alon g th e Fermi line. The shape of the Fermi line 
is shown in ref.ll5|l . 7h os t = ImE^(— i5) for the host system. 



As a result, the structure of the "hot spot" and the "cold 
spot", which are located around (7r/2,7r/2) and (tt/,0) 
respectively, is not smeared out by strong non-magnetic 
impurities. This highly nontrivial result is brought by 
the k-dependence of <5£k- This finding strongly suggests 
that the enhancement of the Hall coefficient near the AF- 
QCP, which is brought by the strong back-flow (current 
vertex correction) around the cold spot [H , does not 
decrease due to the strong impurities. In fact, the Hall 
coefficient for under-doped YBCO slightly increases by 
the doping of non-magnetic impurities |60|. 

On the other hand, in the case of weak impurities 
where the Born approximation is reliable, 7 k mp should 
be almost k-independent. Thus, the structure of the 
"hot/cold spots" will be smeared out by (large num- 
bers of) weak non-magnetic impurities. By this reason, 
the enlarged Hall coefficient near the AF-QCP, which is 
brought by the back-flow around the cold spot, is reduced 
by weak non- magnetic impurities |6lj. This theoretical 
result would be able to explain the reduction of i?n in 
CeCoIns at very low temperatures [4(J, |41j . 

Another important finding is that 7™ p given by the 
GV 7 -method becomes larger than 7hn P ; due to the non- 
local scattering (non s-wave scattering) given by <5£. Fig- 
(77) ure IT31 shows 7™ p (0) for 7i imp = 1 given by the GV 1 - 
mcthod, at T = 0.02 for LSCO, YBCO and NCCO. 
In both hole and electron-doped systems, the hot /cold 
spot structure in the host system remains even in the 
presence of impurities. The absolute value of J^ p (0) 
increases drastically when the system is close to the AF- 
QCP. as shown in^| This finding gives the explanation 
for the huge residual resistivity in metals near the AF- 
QCP, which has been a long-standing problem in strongly 
correlated electron systems. In contrast, 7™ P given by 
the GV°-method is comparable to 7i° np , because the en- 
hancement of x s around the impurity is not taken into 
account. Also, 7™ p by the Gl^-method is much smaller 
than 7 mip . This result will be an artifact of the GV- 
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FIG. 14: (Color online) j™ PW (I = 1 ~ 4) per n imp at 
T = 0.02, along the Fermi line. They are given by the GV 1 - 
method. They represent contributions by eas, 17H - 174L re- 
spectively. 



as U is increased. In the same way, <y™ p( - 4 - > takes a large 
value for NCCO at lower temperatures. This enhance- 
ment of ^™ p ( 4 ' at lower temperatures gives rise to the in- 
sulating behavior of the resistivity, as we will show later. 



FIG. 13: (Color online) 7^(0) per n imp at T — 0.02, along 
the Fermi line. They are given by the GV^-method. Both 
7^ mp and 7h os t show similar k-dependence, which means that 
the hot /cold spot structure is maintained against the doping 
of strong nonmagnetic impurities. Moreover, 7™ p is much 
larger than 7 ; mp when AF fluctuations are strong. 



method, as discussed in Wl 

Here, we examine the origin of the enhancement of 
7™ p in more detail. Figure d shows 7 ™ p(1) (I = 1 ~ 4) 
for LSCO and YBCO, which represent contributions by 
eas. l(7T|) - (|74l . respectively. They are also expressed in 

Fig. EU Note that 7 ™ p = J^Li ^ m ■ In both YBCO 
and LSCO, ^/™ p ^ 2 ^ and ^™ p ( 4 ) give main contributions: 
The latter is dominant for YBCO, whereas 'y™ p ( 2 - ) is com- 
parable to 7 ™ p(4) around k = (7r/2,7r/2) for LSCO. In 
both systems, ^™ p ( 4 ) grows drastically below T — 0.02 



B. Resistivity 

Here, we calculate the resistivity in nearly AF metals 
in the presence of strong impurities. Hereafter, we take 
account of the impurity effect only up to 0(n; mp ). In 
other words, we neglect the interference effect (e.g., the 
weak localization effect) which is given by the higher or- 
der terms with respect to nj mp . In fact, many anoma- 
lous impurity effects of HTSC in which we are inter- 
ested are of the order of 0(n lmp ). For example, the 
residual resistivity in HTSC is proportional to the im- 
purity concentration for n lmp <, 4% 36] . Surprisingly, 
the residual resistivity per impurity increases drastically 
as the system approaches the half- filling. The relation 
Ap ~ {Ah/e 2 )n lmv /5 holds in the under-doped region, 
which is n/S times larger than the residual resistivity in 
2D electron gas. (6 = |1 — n| is the carrier doping con- 
centration.) Hereafter, we will explain this experimental 
fact based on the idea that the effective cross section of 
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FIG. 15: Several diagrams for the conductivity in the presence 
of impurities. The residual resistivity Ap given by type (a) 
and (c) diagrams is proportional to 0(ni mp ). On the other 
hand, Ap given by type (b) diagrams is 0(nf mp ). In the 
present study, we take only type (a) diagrams into account. 



an impurity is enlarged due to many body effect near the 
AF-QCP. 

The conductivity is given by the two-particle Green 
function; we show some diagrams in Fig. If 51 Here, the 
cross represents the impurity potential, and the filled cir- 
cle is the three point vertex due to the electron-electron 
correlation. In Fig. the type (a) diagrams give the 
correction of the order of 0(ni mp ), whereas the type (b) 
diagrams are 0(nf ) because they contain cross terms 
between different impurities [58|. Thus, we drop all the 
diagrams which contain cross terms. The type (c) dia- 
grams, which contain current vertex corrections (CVC's) 
due to impurities, which are necessary to cancel the ef- 
fect of forward scattering on the conductivity. However, 
we drop all the diagrams with CVC's for the simplicity of 
the calculations, which would be allowed for a qualitative 
discussion. We also drop the CVC due to the electron- 
electron correlation because it gives merely a small cor- 
rection to the conductivity, as shown in ref. [l5j . 

As a result, by neglecting the CVC's, the conductivity 
without and with impurities up to 0(nj mp ), op and <7i mp , 
are given by the following equations [SH l58l l59| : 



00 



= e2 E / | \Gi(e)\ 2 vUe), (78) 

k 

= [-(-¥) \G k (e)\ 2 vl x (e), (79) 



dc 



where v kx (e) = dc k /dk x +ReEk(e). G k (e) = (e + ^-e k - 
Sk(e)) -1 is the Green function obtained by the FLEX 
approximation without impurity. The averaged Green 
function in the presence of impurities, G k , is given by 



G k (e - iS) = ({G° k (e - iS)}- 1 - i 7k mp (e)) 



(80) 



Figure 1161 shows the temperature dependences of the 
resistivities for LSCO, YBCO and NCCO obtained by the 
Gy 7 -method; p imp = l/cr imp for n imp = 0.02. p = 1/<tq 
is the resistivity without impurity. We see that the 
"residual resistivity" at finite T, which is the incre- 
ment of resistivity due to impurity Ap = pj mp — po, is 
approximately constant for a wide range of T, except 
for the abrupt increase at lower temperatures in YBCO 
and NCCO. Ap grows as we enlarge U (YBCO, LSCO) 
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FIG. 16: (Color online) p- lmp (p° mp ) represents the resistivity 
for rij m p = 0.02 given by the GV 1 -method (local scattering 
approximation), po is the resistivity for the host system. At 
higher temperatures, the residual resistivity pi mp — po is almost 
T- independent, and its value increases near the AF-QCP. At 
lower temperatures, the insulating behavior emerges in the 
close vicinity of the AF-QCP. T = O.f and p= 1 correspond 
to 400K and 250pf2-cm for HTSC's, respectively. 
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1 < 7° 



or change the filling number n towards the half-filling 
(NCCO). This result comes from the fact that the scat- 
tering cross section of an impurity is enlarged by the non- 
local modulation of the self-energy, ST*. In other words, 
non s-wave (elastic and inelastic) scatterings caused by 
(5£ gives an anomalously large residual resistivity near 
the AF-QCP. 

Figure El also shows that pf mp derived from eq. i(7$|) 

by replacing 7 k mp (e) in eq. (JHOJ) with 7 i t J np (e) which is 
momentum independent. We call it the "local scatter- 
ing approximation" because non-s-wave impurity scat- 
tering processes caused by the nonlocal effect (<5I!) are 
dropped. In contrast to Ap, the doping dependences of 
Ap° = pf mp — po is much moderate. Moreover, Ap° 
slightly decreases at low temperatures as 7k becomes 

anisotropic: In fact, ( „ — )~ 1 — ( — ) 

7k + 7iV V 

when 7k(^ 7hn P ) i s ver y anisotropic. 

In each compound (YBCO, LSCO, NCCO), the av- 
erage spacing between Cu02-layers is about 6A. Using 
the relation h/e 2 = 26kf2, p = 1 in the present calcu- 
lation corresponds to 250/if2-cm. As shown in fig. El 
residual resistivities for nj mp = 0.02 at T = 0.05 for 
LSCO (U = 5), YBCO (U = 8) and NCCO (n = 1.13) 
are 0.7, 1.0 and 0.8, respectively. They correspond to 
175 ~ 250pft-cm. These obtained values are close to the 
experimental residual resistivity in under-doped HTSC's; 
Ap~ (4h/e 2 )n imp /\l-n\ ^. 

Another important findings is the insulating behavior 
(dp imp /dT < 0) in YBCO and NCCO at lower temper- 
atures. We also checked that a similar upturn is also 
observed in LSCO for U = 6. This insulating behav- 
ior of pi mp is caused by the steep increment of 7 k mp at 

lower temperatures, which is mainly given by <y™ p( ' 4 ' ) oc 
lnu)£(0). Therefore, the physical origin of this phe- 
nomenon is the strong inelastic scattering around the im- 
purity. The obtained insulating behavior will be univer- 
sal for systems with strong nonmagnetic impurities in the 
close vicinity of the AF-QCP. Actually, the upturn of p 
is widely observed in under-doped HTSC's, both in hole- 
doped and electron-doped compounds |3al43ll44l4al46]| . 
The origin is not the weak localization as we explained in 
|U Based on the present study, we expect that residual 
disorder will be the origin of the insulating behavior in 
HTSC's. 

In Fig. El we show that U dependences of Ap and Ap° 
(per impurity) at T — 0.02 given by the Gl^-method, 
both for LSCO and YBCO. Ap° increases as U is raised, 
inversely proportional to the DOS given by the FLEX ap- 
proximation, -/Vhost(O). We stress that Ap increases dras- 
tically as U is raised, much faster than Ap° does. This 
prominent increment is derived from the strong inelastic 
scattering around the impurity, Imi5I](0). For compar- 
ison, we also show results given by the GV and GV°- 
methods in Fig. El Ap[ GV0 ] is close to Ap° because the 
enhancement of AF fluctuations around the impurity are 
not taken into account. Even worse, Ap[Gv] decreases 
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FIG. 17: (Color online) Ap = p in 



po per rib 



at T = 



0.02. Ap given by the GV 7 -method grows drastically as U 
increases, whereas Ap by the GV^-method gradually drops. 
Here, "by 7i° mp " represents the residual resistivity by the local 
approximation, which slightly increases with U in proportion 
to 7i ° mp oc l/AT host (0). 



for U > 3 in LSCO. This result should be an artifact 
of the GF-method, because it fails to reproduce the en- 
hancement of spin fluctuations as explained in previous 
sections. 

The prominent (7-dependence of Ap given by GV 1 - 
method would corresponds to the pressure dependence of 
Ap observed in K-(BEDT-TTF) 4 Hg 2 . 89 Br 8 , which stays 
near the AF-QCP at ambient pressure 01 ■ -"- n this com- 
pound, Ap decreases to one tenth of its original value 
by applying the pressure. Considering that the applied 
pressure makes U/Wb&nd small, this experimental result 
is consistent with Fig. El We also note that the residual 
resistivity in heavy fermion systems near the magnetic 
instabilities frequently show prominent pressure depen- 
dence; Ap takes a maximum value at the AF-QCP, and 
it decreases quickly as the system goes away from the 
AF-QCP. H3, HE E3, EH This experimental fact is 
well explained by Fig. 1171 

Figure El shows the filling dependence of Ap given by 
the GT^-method at T = 0.02, which is above the upturn 
temperature of p for YBCO as recognized in Fig. 1101 In 
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FIG. 18: (Color online) Filling-dependence of Ap/ni mp given 
by the GV 7 -method at T — 0.02. Origin of the huge Ap in 
under-doped region is the short r around impurities induced 
by strong AF fluctuations. 



both LSCO and YBCO, Ap increases drastically as n ap- 
proaches unity, far beyond the s-wave unitary scattering 
value (~ 4/n). The obtained result is consistent with ex- 
perimental observations in HTSC's, where Ap ~ Ap° in 
the over-doped region, whereas Ap 3> Ap° in the under- 
doped region [36|. 

In the present calculation, we dropped all the CVC's 
for simplicity. As shown in ref. |l5|. po without impuri- 
ties is slightly enlarged by the CVC; effect of the CVC 
for po is small because the origin of the resistivity is the 
large angle scattering due to AF fluctuations (q w (ir, it)). 
[In general, CVC's for the resistivity are important when 
small angle scatterings are dominant.] In the present 

calculation, the origin of the huge Ap oc <y™ p< ' 4 ' ) is also 
the AF fluctuations induced around the impurity. There- 
fore, we expect that obtained Ap in this section is qual- 
itatively reliable. On the other hand.CVC's play quite 
important roles on i?H, Ap/ p and v j5J]: It is an impor- 
tant future issue to study these transport coefficients in 
the presence of impurities by taking CVC's into account. 
Finally, we note that the obtained insulating behavior 
of Pimp might be under- estimated because the area of re- 
gion A in the present numerical study (M = 6^8) would 
be not enough in the very close vicinity of the AF-QCP 
(£af » 1). 



C. Decrease of the Hole Density (= |1 — n|) 
around the Impurity, and Increase of Ap at T = 

Up to now, we have found that ^™ p ( 4 ' (and 7™ p ^ 2 ' ) ) 
give the main contribution to the huge residual resistivity 
Api m p at finite temperatures. However, 7™ P ^(0) = 
at zero temperature because <5£; jm (0) becomes a real 
function at T = 0. Therefore, it is highly nontrivial 
to predict the value of Ap; mp at T = 0. On the other 



hand, 7^ n , 7™ p ^ and <y™ P ^ could give an enlarged 
Apimp 3> Ap? mp at zero temperatures, because they are 
finite even at T = 0. 

Hereafter, we discuss the residual resistivity at T = 0. 
For simplicity, we assume an two-dimensional isotropic 
system with the dispersion ek — k 2 /2m. At T = 0, 7™ p 
is given by 

imp T rr 

7k = "impImT^k 

2 

= "imp— V"sin 2 (5i, (81) 
i 

where I represents the angular momentum with respect 
to the impurity potential (I = 0, ±1,±2, • • •), and Si is 
the phase shift for channel I. If we drop the CVC, the 
resistivity at T = is given by 

Pimp = ^?E Sin2 ^- ( 82 ) 

I 

Here, sin 2 Si is replaced with sin 2 Si — cos (Si — 
Si ) sin Si sin Si if the CVC due to impurity is taken 
into account. On the other hand, the number of local- 
ized electrons around the impurity, Antot , in an isotropic 
2D system is given by j6^ 

An to t = ~YSi. (83) 

I 

Thus, the residual resistivity at T = will grows as the 
number of electrons in bound states increases. 

YBCO 



T 


0.01 


0.015 


0.02 


0.05 


100 x An(l,0) 


6.15 


1.50 


1.06 


0.60 


Antot 


2.86 


0.334 


0.190 


0.052 


ast 


0.9968 


0.9950 


0.9932 


0.9807 



NCCO 



T 


0.02 


0.025 


0.03 


0.06 


100 x An(l,0) 


-3.84 


-3.12 


-2.78 


-2.13 


Antot 


-1.52 


-0.897 


-0.639 


-0.254 


ast 


0.9965 


0.9958 


0.9948 


0.9839 



TABLE I: Temperature dependences of An(l,0), Antot and 
the Stoner factor ast = maxqt/Il^O). 

Table [I] represent An(l,0) = n(l,0) — n and Antot = 
J2 r An(r) for both YBCO (n = 0.9, U = 8) and NCCO 
(n = 1.13, U = 5.5). We see that the number of elec- 
tron at r = (1,0) approaches the half filling (n = 1) 
at low temperatures for both YBCO and NCCO. More- 
over, An to t increases (decreases) prominently in YBCO 
(NCCO) as T decreases. This result suggests that the 
residual resistivities at T = will be larger than pp mp in 
the vicinity of the AF-QCP. 
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Finally, we discuss why n(l,0) approaches unity 
and An tot increases (decreases) in hole-doped (electron- 
doped) systems. In the FLEX approximation, the ther- 
modynamic potential Q in a uniform system is given by 

n = -T^Tr [SG + ln(-[G ]- 1 +53)] 



+T^Tr 
i,l 

+UU a + U 2 [U ] 2 ] 



|ln(l - UTl°) + -ln(l + ?/n n ) 



(84) 



As we explained, the AF fluctuations are enhanced 
around the impurity. This effect would be expressed by 
increasing II (q, u>i) in eq. ijSijl by An (> Q), or U by 
AC/ (> 0). According to eq. JH3}, 

W ra -lrX;n («Lw«)(i-tm (q lWI ))- 1 > (85) 



in the case of 1 - £/n°(Q,0) < 1. In deriving (pjj) . 
we used the fact that the implicit derivative through S 
vanishes because of the stationary condition <557/<5£ = 
in the conserving approximation [(J, l64j . According to 
eq.lJSSJ, we obtain 



On 
dU 



dU \dfi 



c>II (q,u; ; ) 



2 ^(l-C/IT°(q, W; )) 2 dp 



(86) 



As a result, the electron density n around the impurity, 
where AU > is satisfied as mentioned above, increase 
(decreases) when dU(Q,0)/dp is positive (negative). 
Therefore, we conclude that An tot will increase (de- 
crease) in hole-doped (electron-doped) systems, as rec- 
ognized by numerical results given by the Gl^-method. 



VII. DISCUSSIONS 



A. Summary of the Present Work and Future 
Problems 



The present study reveals that a single impurity 
strongly affects the electronic states in a wide area around 
the impurity in the vicinity of the AF-QCP. For this pur- 
pose, we developed the GV^-FLEX method, which is a 
powerful method to study the impurity effect in strongly 
correlated systems. The GV 1 method is much superior 
to the GV, which is a fully self-consistent FLEX approx- 
imation. Using the GV 1 method, characteristic impurity 
effects in under-doped HTSC's are well explained in a 
unified way, without introducing any exotic mechanisms 
assuming the breakdown of the Fermi liquid state. The 
main numerical results are shown in Fig. [S] (local DOS 



around the impurity site), Figs. 16181 flocal and staggered 
susceptibilities), and Figs. 1161181 (resistivity in the pres- 
ence of impurities). Qualitatively, these obtained numer- 
ical results are very similar for YBCO, LSCO and NCCO. 
Therefore, novel impurity effects in nearly AF metals re- 
vealed by the present work would be be universal. 

Based on the GV 1 method, we found that both lo- 
cal and staggered susceptibilities are prominently en- 
hanced around the impurity site, as shown in Figs. EJand 
Especially, a nonmagnetic impurity causes a Curie- 
like spin susceptibility, p 2 s /ST. The GU / -method gives 
p cH ~ 0.74^ B for LSCO (n = 0.9, U = 5), shown in Figs. 
03 Note that /Lt eff - lfi B in YBa 2 Cu 3 6 . 66 (T c » 60K). 
We also found that the quasiparticle damping rate takes 
a huge value around the impurity, owing to the enhanced 
AF fluctuations. By this reason, the local DOS at the 
Fermi level is strongly suppressed around the impurity 
site, which forms a so-called "Swiss cheese structure" as 
shown in Fig. ^3 Its radius is about the AF correlation 
length of the host, £af, which is about 3 ~ 4a (a be- 
ing the lattice spacing) in slightly under-doped HTSC's 
at T = 0.02. We guess that Swiss cheese holes stay in 
a normal state even below T c , because of the extremely 
short quasiparticle lifetime there. In fact, the residual 
specific heat (T « r c ) induced by an impurity becomes 
very large in under-doped systems |50j . This experimen- 
tal fact will be explained in our future study based on 
the GV^-mcthod 

Near the AF-QCP, the short quasiparticle lifetime in- 
side the Swiss cheese hole gives rise to a huge residual 
resistivity Ap, as shown in Figs. 1 1 61 1 81 In the under- 
doped region, Ap grows far beyond the s-wave unitary 
scattering limit ~ (4Ti/e 2 )ni mp /n. We find that Ap is 
almost T-independent for a wide range of temperature, 
and it increases drastically as the system approaches the 
AF-QCP. This result is consistent with experiments for 
HTSC's. The obtained value of Ap, 175 ~ 250/xfi-cm for 



= 0.02, are recognized in under-doped HTSC's \3t 



Furthermore, in the close vicinity of the AF-QCP, the 
resistivity given by the Gl^-method shows the "Kondo- 
like" insulating behavior (dp/dT < 0) under the pres- 
ence of nonmagnetic impurities with low concentration 
(~ 2%). This surprising result would explain the "upturn 
of resistivity" which is frequently observed in under-dope 
HTSC's, by assuming the existence of residual disorders. 
The mechanism of this insulating behavior had been a 
long-standing unsolved issue in under-doped HTSC's. 
Different from a conventional single-channel Kondo ef- 
fect, the residual resistivity given by the GU 7 -method 
grows far beyond (Ah/e 2 )ni mp /n. 

We stress that that 7™ P given by the GT^-method 
has strong k-dependence, so the structure of "hot/cold 
spots" is maintained against impurity doping. This re- 
sult could be examined by the ARPES measurements. 
We also comment that the CVC's due to spin fluctua- 
tions cause a finite "residual resistivity", if we define it 
as the extrapolated value at T = 0, even in the absence 
of impurities . We have to take this fact into account 
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FIG. 19: Swiss cheese structure induced by nonmagnetic im- 
purities (nj mp = 0.021) denoted by open circles, which replace 
Cu sites randomly. In the shadow area, both local and stag- 
gered spin susceptibilities are enhanced, and the quasiparticle 
lifetime is depressed. The radius of the Swiss cheese hole is 
~ £af- 



in analysing experimental data. 

In the FLEX approximation, the AF-order (in the 
RPA) is suppressed by thermal and quantum fluctua- 
tions, which is expressed by S°. In the Gl^-method, 
the reduction of fluctuations around the impurity site 
gives rise to the enhancement of susceptibility. However, 
this mechanism is absent in the RPA. Therefore, the en- 
hancement of susceptibility is tiny within the RPA. In 
addition, in the G^ 7 -method, the spin (charge) suscepti- 
bility x Is ( c ) contains the self-energy correction by U's and 
that by I's are treated on the same footing, whereas the 
cross term <5£ is dropped because it will be cancelled out 
by VC's in large part. On the other hand, ^ RPA given 
by the RPA contains only self-energy correction by 7's 
HE HH • This equal footing treatment of the corre- 
lation effect and the impurity effect in the G^ 7 -method 
would be the reason for the superiority of this method. 

In future, we will study the transport phenomena by 
taking the current vertex corrections (CVC) accurately, 
in order to solve the impurity effect on various trans- 
port coefficients in HTSC's and in related systems. [6l| . 
We will also study the superconducting state around a 
nonmagnetic impurity in under-doped HTSC's |6lj . to 
explain experimental observations given by STM/STS 
measurements |48l lo lL l53 |. As for HF systems and in 
organic metals, the strength of residual impurity (or dis- 
order) potential would be comparable with bandwidth 
W^band- Therefore, we have to study to what extent the 
obtained results in the present paper (for I = oo) hold 
in the case of I ~ Wband- 



B. Possibility of the Impurity-Induced Magnetic 
Order 

In the GF 7 -method, the self-energy for the host system 
is given by the FLEX approximation, which cannot be 
applicable for heavily under-doped systems near the Mott 
insulating state. However, FLEX approximation gives 
qualitatively reliable results for slightly under- doped re- 
gion (|l-n| ~ 0.1) to over-doped region HHEIlllll. 
To obtain quantitatively reliable results for under-doped 
region, vertex correction (VC) for the self-energy will be 
necessary, as indicated by ref. [65| . The pseudo-gap phe- 
nomena under T* are well reproduced by the FLEX+T- 
matrix approximation, where self-energy correction due 
to strong superconducting (SC) fluctuations are taken 
into account 0, H El 13. By taking the CVC due 
to SC fluctuations into consideration, anomalous trans- 
port phenomena in HTSC's under T* (e.g., prominent 
enhancement of the Nernst coefficient) are well under- 
stood in a unified way p| . 

One of the merits of the FLEX approximation is that 
the Mermin- Wagner theorem with respect to the mag- 
netic instability is satisfied. In fact, previous numerical 
studies based on the FLEX report that no SDW-order 
emerges in two dimensional systems at finite T. In Ap- 
pendix A, we offer an strong analytical evidence that the 
FLEX approximation satisfies the Mermin- Wagner theo- 
rem. Also, the GV^-method does not predict any SDW- 
order in two dimensional systems with single impurities, 
at least for model parameters studied in the present pa- 
per. However, when the concentration of impurities is fi- 
nite, Swiss cheese holes will overlap when £af exceeds the 
mean distance between impurities (I). Therefore, a SDW- 
order or a spin-glass order would happen when £af > I- 
In fact, in Zn-doped YBCO, the freezing of local mo- 
ments due to Zn was observed by /u-SR study at very low 
temperatures |66| . Note that impurity induced AF order 
occurs in two-leg ladder Heisenberg models ■ 



C. Comments on Related Theoretical Works 

Here, we discuss the impurity effect on the electronic 
states of the "host system" . We did not study this effect 
because it is higher order effect with respect to n; mp . In 
Zn-doped YBa 2 Cu 4 8 , Ihot et al. measured 63 Cu 1/TjT 
in the host system (i.e., away from the Zn sites) [6c| . 
and found that the host AF fluctuations above T* are 
reduced by a few percent doping of Zn. They suggest 
that the localization effect would be the origin of this 
reduction. 

Based on the FLEX approximation (GV-method in the 
present paper), authors of ref. 69] studied the impu- 
rity effect on the electronic state by neglecting the k- 
dependence of the self-energy. They reported that the 
AF fluctuation in the host system is depressed by impu- 
rities, which seems to be consistent with ref. [6^. As we 
have shown however, the Gl^-method fails to reproduce 
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correct electronic states around the impurity: In fact, 
both the Curie like A% s as well as the huge Ap are satis- 
factorily reproduced only by the GV^-method. We note 
that p{T) in HTSC shows a nearly parallel shift by the 
impurity doping, which means that the inelastic scatter- 
ing in the host system is independent of impurities. This 
fact would suggest that the AF fluctuations in the host 
system is unchanged, in contrast to the NMR result [68| . 
It would be an important issue to understand these two 
experimental facts consistently. 

Next, we comment that Miyake et al. have intensively 
studied the impurity effects near QCP • They found 
that Ap due to weak impurities (Born scattering) is en- 
larged when the charge fluctuations are developed. Their 
analysis corresponds to GV°-method since the suscepti- 
bility is assumed to be independent of impurities. In 
contrast, the present work based on Gl^-method shows 
that the Ap due to strong impurities increases near the 
AF-QCP, which originates from the enhancement of x Is 
around the impurity. We note that Xflex = dn/dfi is 
less than half of its non-interacting value in the present 
range of parameters. In analysing experiments, we have 
to consider carefully which kind of criticality may occur 
in the compound under consideration. 

Finally, we explain previous studies on disordered Hub- 
bard model based on the dynamical-mean-field-theory 
(DMFT) in the d = oo limit [H El 113 , and we make 
comparison with the present study. In the DMFT, the 
local-moment formation is found when strong hopping 
disorders (offdiagonal disorders) exist Also, it is 

found that the Neel temperature increases due to weak 
onsite disorders (diagonal disorders); I < U Jz^. How- 
ever, a local-moment formation outsside of the impu- 
rity and a huge residual resistivity, which are realized 
in HTSC's, cannot be explained by the DMFT. For this 
purpose, a nonlocal modulation of the self-energy around 
the impurity have to be taken into account, which is pos- 
sible in the GF 7 -method. 
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APPENDIX A: MERMIN- WAGNER THEOREM 
IN 2D ELECTRON SYSTEMS 

The Mermin- Wagner (M-W) theorem states that any 
magnetic instabilities are absent at finite T in 2D sys- 
tems. Actually, the SCR theory satisfies the M-W theo- 
rem 0. As for the FLEX approximation, however, the 
M-W theorem had been recognized only by numerical 
studies. In this appendix, we present a strong analytical 
evidence that the FLEX-type self-consistent spin fluctu- 
ation theory satisfies the M-W theorem. 

Here we introduce the phcnomenological expression for 
the dynamical spin susceptibility as 

XqH= i + Ci F (q-Q) 2 --M f ' (M) 

where £af is the AF correlation length. Q is one of 
the nesting vectors which minimize | q — Q | . Apparently, 
Xq(0) = X-q(0)- Both xq an d w^ 1 are proportional to 
£af m the FLEX approximation. 

When the system is very close to the AF phase at finite 
temperatures where lo s { 3> T is satisfied, the self-energy 
within the scheme of the FLEX approximation is given 
by 

E k (iw„) ps T Gk+q(w n )Xq(0) 
q 

w G k+Q {iLu n )A (A2) 

Q7/2 

A = T— £ Xq (0) (A3) 
q 

where the static approximation is applied, which will offer 
the upper limit of Tjy. A diverges as T — > Tn in propor- 
tion to T In £af • Especially, A = oo at T = T N (> 0). 
On the other hand, A is finite even at T = Tjy in 3D 
systems. According to ca. (|A2|) . 



iU„ — £k+Q — Ek+Q^n) 



lUJ n - £ k +Q - iaJ „_ ek _E k (i w „) 

Equation l|A4p can be solved analytically. Considering 
Ek(ct>) = when A = 0, the self-energy and the Green 
function for real frequencies are given by 



E k (w) = -^-ek)-sgn(^-e k )^/>-e k ) 2 -4A - £k (A5) 



G k ( W ) = i-(oj- e k ) + sgn(c - e^-Jiu, - 6 k ) 2 - AA - ^ (A6) 

\ 2 1 y w - £ k+Q j 

One can check that u> ■ ImSk(w) < 0. This is not a Fermi liquid because the renormalization factor z is zero (owing 



-l 
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to the static approximation). Hereafter, we assume ek = 2t(cosk x + cosk y ) at half filling (n = 1), that is, both 
the perfect nesting and the particle-hole symmetry exist. Apparently, Q = (tt,tt), £k+Q = — e ki and fJ> = 0. The 
irreducible susceptibility at q = Q and to = is given by 

n Q(°) = -£ / ^th^M^ +Q HG^( W )} (A7) 

k ^ 

According to eq. I|A6(1 . 

G k+Q (u;)G k (u;) = Q(^ 2 - e k ) + isgn(^ 2 - e 2 ) V^H + Jsgn(^ 2 - e 2 )|^ 2 - e 2 - 4A|) (A8) 
5k M = (u; 2 -e k )(u; 2 -e k -4A) (A9) 
<7 k (u>) is negative when |ek| < \oj\ < a/c^. + 4A. Apparently, the integrand in ea. (|A7fl is finite only when g k (w) < 0. 

- Im{G k+Q (a;)G k (w)} = -Im{( A + ^y/-gk(w) T 1 } < for 5k (w) < (A10) 

= for g k (u) > (All) 



According to ea. (|A7(l . in the case of \f~A Wband, 

n Q (0) ~ OiA- 1 ' 2 ) (A12) 

Note that ImGk(w) is non-zero when ffk(w) < 0, and 
one can check that dujlmG^(uj) = —n for any 
A by MATHEMATICA°° Considering that xq(0) = 
n Q (0 )/(l - [/IIq(Q)) in the FLEX approximation, eq. 
(|A12|) means that Xq(0) approaches to zero as T — ► T)v 
(because A — > oo in the case of TV > 0) in 2D sys- 
tems. Equation 1JA12I) suggests that £af oc e 1 / 7 when 
the ground state is a ordered state. 

As a result, TV cannot take a finite value in the two 
dimensional model with perfect nesting at half filling, 
which is the most likely to cause the magnetic instabil- 
ity. Therefore, the present analysis gives a compelling 



evidence the the FLEX approximation satisfies the M-W 
theorem for general 2D systems. 

Here we rewrite A as A' + A", where A' = 
3 jj2 lq-QI<9c 

T ^ Xq(0), and q c is a cutoff momentum, 

q 

Then, A" is a smooth function around the AF-QCP. 
When T > + q% , A' will be the number of free bosons 
(magnons) with energy e q = q 2 within the radius of q c 
at the chemical potential /i = — £^f- Thus, a magnetic 
instability (£af — > oo) corresponds to the Bose-Einstein 
condensation of magnons (fi — 0). Therefore, a meaning- 
ful correspondence between the M-W theorem and the 
absence (presence) of the Bose-Einstein condensation in 
2D (3D) systems is recognized. 
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